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A new approach to multi-jet calculations in hadron collisions 

M. Moretti a 

a Department of Physics, University of Ferrara and INFN sezione di Ferrara 

The ALPHA algorithm to evaluate the exact, tree- level matrix elements is reviewed in the context of multi- 
parton processes in QCD. The algorithm is suited for the authomatic calculation of tree- level scattering amplitudes 
and allows for straightforward inclusion of mass effects. It's CPU cost grows like K n (K constant, n number of 
external particles) as opposed to the factorial growth of Feynman graphs. I also discuss the summation over colour 
configurations which is designed to allow the construction of parton-level event generators suitable to interfacing 
with a parton-shower evolution including the effects of colour-coherence. Explicit results for the total rates and 
differential distributions of processes with 8 final-state partons are given. 



1. Introduction 

The ability to evaluate production rates for 
multi-jet final states will be fundamental at the 
LHC to study a large class of processes, within 
and beyond the SM. A necessary feature of any 
multi-jet calculation is the possibility to prop- 
erly evolve the purely partonic final state, for 
which exact fixed-order perturbative calculations 
can be performed, into the observable hadronic 
final state. This evolution is best performed us- 
ing shower Monte Carlo calculations. The ac- 
curate description of color-coherence effects, fur- 
thermore, requires a careful bookkeeping of the 
contribution to the matrix elements of all possible 
color configurations. The goal of the algorithm [ 
|l| described in this talk is to allow the effective 
calculation of multi-parton matrix elements, al- 
lowing the separation, to the leading order in 
1/N% (N c = 3 being the number of colors), of the 
independent color configurations. This technique 
allows an unweighting of the color configurations, 
and allows the merging of the parton level calcu- 
lation with the HERWIG Monte Carlo. 

The key element of the strategy is the use of 
the algorithm ALPHA, introduced in Ref. [ |2| for 
the evaluation of arbitrary multi-parton matrix 
elements. The algorithm is built up from an it- 
eration of matrix multiplication and is intended 
for the authomatic calculation of tree-level am- 
plitudes. It has a complexity growing like a 
power in the number of particles, compared to 



the factorial-like growth that one expects from 
naive diagram counting. 

The ALPHA algorithm will be reviewed in section 
(|^), and its application to multi-jets physics in 
section (||). 

2. The ALPHA algorithm 

In reference [ 0], a new approach to the com- 
putation of tree level scattering amplitudes was 
introduced. This approach, based on the numer- 
ical Legendre transform of the effective action, is 
particularly useful for the authomatic calculation 
of multi-particle processes. This technique was 
implemented in a FORTRAN code [ || which has 
been succesfully used to study several intricated 
electroweak processes [ 0, [| of interest both at 
LEP [ [| and at the NLC Among the key fea- 
tures of the algorithm are: 

• Suitable (easy to implement) for the au- 
thomatic calculation of the scattering am- 
plitudes. 

• Inclusion of mass effects is straightforward 

• The CPU cost grows like K n (K constant, n 
number of external particles) as opposed to 
the factorial growth of the number of Feyn- 
man graphs. 

We review in this Section the ALPHA algorithm; 
the interested reader can find a more detailed dis- 
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cussion in the original paper [ [l], || , which includes 
an explicit analytic example for the A<fi 3 theory. 

Let r be the one-particle-irreducible generator 
of the Green functions for a given theory. Then 
the computation of the S-matrix requires the eval- 
uation of the Legendre transform, Z, of F: 

z(j a ) = -r(0 Q ) + j a (x)0 a (x) 

where (f> a are the classical fields defined as the 
solutions of 



J a = 



ST 



and the J a play the role of classical sources. 

For concretness we will develop in some detail 
the case of the scattering amplitue for n exter- 
nal gluons. At tree- level T coincide with the La- 
grangian. In momentum space 



-l/2( Pfl A«- Pl/ A«) 2 
+9f ab c{p»Al-p v Al)AlAl 
-(B;j 2 - 2gf ahc Bl v A\Al + J£A% 



where B^ v is an auxiliary field and it is intro- 
duced in order to deal with trilinear interactions 
only. The sources J° is a standard source term: 
J(p) = Y^j=i e< ^,HP~ Pj)i i- e - ^ contains the rele- 
vant excitation for the external particles. Notice 
that because of this choice only a finite number of 
momenta (linear combinations of pj) enters into 
the problem which is reduced to a problem with 
a finite number of degrees of freedom. 

The fields A(p) are then found as solutions of 
the equation of motion. Let us stick to Feynman 
gauge for definitness. 



y abc [2(k-p)-A b (q)A^(k) 



~2q,A b (q)-A^k) 
-Bl( q )A c Jk)]+^J,(p) 
B"(p) = -gf ab cA b Aq)Al{k) 



(1) 



where momentum conservation is understood. 

The equations of motion (|l|) are solved itera- 
tively (expansion in g) and this implies that the 



problem is solved with a loop of matrices multi- 
plication, more suitable for numerical implemen- 
tation than the standard approach. The initial- 
ization steps are 



= o 



(2) 



and the subsequent steps 



A;(Pj+Pk) 

gfa 



be 



(Pj +Pk) 2 



[-2{p 3 +2p k )-A b {p j )A c tl {p k ) 

+2( Pk + 2p 3 )- A b {p k )A c )1 {p 3 ) 
-2fe+p fe ) AI A b (p fc ).A c (^) 



B%{Pj+Pk) = -gfa b cA b (p k )A^p 3 ) 



(3) 



The step (g|) is then iterated until we have con- 
structed A and B fields with up to n momenta. 
An important remark is in order here. Per- 
forming the iteration (^) we drop terms contain- 
ing e a (pj)e a (pj) i.e. twice the same external mo- 
menta. In fact, although legitimate, they don't 
contribute to the amplitude. 
Finally the amplitude is 



(4) 



A(p 3l ,. --,Pj n ) = Cym(A) 



where the A fields are obtained in eqns. 
Notice that the prescription to drop terms con- 
taining at least twice the same external momenta 
(~ e a (pj)e a (pj)) is still kept and it is because of 
this prescription toghether with the choice (||) for 
the initial step that neither truncation nor func- 
tional derivation is actually required. 

A remark is in order here. The algorithm 
sketched above has an important feature: it pro- 
vides a very compact way of storing the relevant 
information. Indeed the number of contraction 
A a (pj 1 +. . -+Pj m ) (see eq. ||) which is needed is of 
the order 2", each of them requiring roughly the 

1 Using the equation of motion it is possible to halve the 
required number of iteration step as well as to reduce the 
number of contributions entering the final expression of 
the amplitude. 
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same CPU time to be computed. Therefore the 
CPU cost of the algorithm grows like a constant to 
the power n, n being the number of external par- 
ticles, as opposed to the factorial growth of the 
number of Feynman graphs. To hint how does 
this work let us develop in more detail the case of 
the amplitude for five external gluons. For sim- 
plicity we neglect the four gluons coupling. The 
A fields are given by 

A a M) = *>i) 3 = 1,..., 5 

, g {°% ["2(Pi + 2p fe ) • A»( Pj )A^p k ) 

+2(p k + 2 Pj ).A b ( Pk )A^( Pj ) 
-2( Pl +p k )^A b (p k )-A c ( Pj )} 
j = 1, . . . ,5 

(5) 

and the ALPHA amplitude by 

5 

A = gJ2 Aa (Pi)° 

3=1 

5 

^ fabc.T a l3 1 A b ( Pk + P i)pA C ( P 

m i P n ) 7 

k^l^m^n^j 

T a f3~y = g a p(pj ~ P k - P l)~/ 

■in + Pn -Pj)fi 
+9/3j(Pk + Pi - Pm ~ Pn)a 

(6) 

notice that the summations over k, I, m, n, as well 
as over lorentz and color indices, are carried over 
before the multiplication by A a (pj). It is this 
feature extensively used in the construction of 
both the A fields and the amplitude which turns 
the factorial growth of the number of Feynman 
graphs to a power law growth of the CPU cost of 
the ALPHA algorithm. Notice that in eqns. (||J6|) 
A(pj + Pk + Pm) and A(pj + Pk + Pm + Pi) are 
not computed (and used). Morover neither the 
kinetic term nor the source one appears in the 
expression (0) of the amplitude. These simplifi- 
cations occur because of the equation of motion 
as discussed in more detail in [ 0, 0| • 



3. Multi-jets processes 

Multi-jet final states play an important role in 
the study of high-energy collisions. They provide 
in fact interesting signatures for several phenom- 
ena, both within the Standard Model (e.g. top- 
pair production), and beyond it (e.g. multi-jet de- 
cays of supersymmetric particles such as gluinos 
and squarks). The accurate determination of the 
properties of these phenomena requires a good 
understanding of the properties of the usually 
large multi-jet QCD backgrounds, which can dis- 
tort the shapes of signal distributions and affect 
the measurement of quantities such as resonances' 
masses. 

There are several reasons for wanting to im- 
prove the tools currently available to perform 
these calculations. 

1. First of all, interesting final states with 
larger jet multiplicities will become avail- 
able with the next generation of colliders 
(LHC and NLC). 

2. Secondly, one would like to be able to com- 
plement the calculation of parton-level ma- 
trix elements with the evaluation of the full 
hadronic structure of the final state. 

The goal of the algorithm [ [jj described in this 
Section is to allow the effective calculation of 
multi-parton matrix elements, allowing the sep- 
aration, to the leading order in 1/N% (N c = 3 
being the number of colors), of the independent 
color configurations. This technique allows an 
unweighting of the color configurations, and al- 
lows the merging of the parton level calculation 
with the HERWIG Monte Carlo. The key ingredient 
of our strategy is the ALPHA algorithm outlined 
in the previous section. It's power-law growing 
in complexity is a necessary feature of any at- 
tempt to evaluate matrix elements for processes 
with large numbers of external particles, since the 
number of Feynman diagrams grows very quickly 
(see table 1) beyond any reasonable value. 

Once the hard scattering matrix element is 
known the prescription [ to correctly generate 
the parton-shower associated to a given event in 
the large- N c limit is the following: 
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Process 


n = 7 


n = 8 


n = 9 


n = 10 


9 9^ng 


559,405 


10,525,900 


224,449,225 


5,348,843,500 


qq^ng 


231,280 


4,016,775 


79,603,720 


1,773,172,275 



Table 1 



Number of Feynman diagrams corresponding to amplitudes with different numbers of quarks and gluons. 



1. Calculate the (n — 1)1 dual amplitudes cor- 
responding to all possible planar color con- 
figurations. This can be done with ALPHA 
using N c large enough [ [l|] . 

2. Extract the most likely color configuration 
for this event on a statistical basis, accord- 
ing to the relative contribution of the single 
configurations to the total event weight 0. 
Since each dual amplitude is gauge invari- 
ant, the choice of color-configurations is also 
a gauge-invariant operation. 

3. Develop the PS out of each initial and final- 
state parton, starting from the selected 
color configuration. This step can be car- 
ried out by feeding the generated event to 
a Monte-Carlo program such as HERWIG, 
which is precisely designed to turn partons 
into jets starting from an assigned color- 
ordered configuration. 

Notice that, if the dual amplitudes are evaluated 
for a specific helicity configuration, HERWIG will 
also include spin-correlation effects in the evolu- 
tion of the parton shower [ ||, ||, @, |[ D . 

As a result, use of the dual-amplitude repre- 
sentation of a multi-gluon amplitude allows to 
accurately describe not only the large-angle cor- 
relations in multi-jet final states, but also the full 
shower evolution of the initial- and final-state par- 
tons with the same accuracy available in HERWIG 
for the description of 2-jet final states. 

In alternative to the above prescription, one 
can use ALPHA to calculate the matrix elements 

2 Defining wt = Aj| 2 for each color flow i, and Wi = 
Zjfc— i ' Wk /^2k— l n w k> •'"th c °l° r structure 

will be selected if Wj—i < £ < Wj, for a random num- 
ber £ uniformly distributed over the interval [0, 1]. 



for external states with assigned colors. Since 
these states are all orthogonal, such an approach 
is particularly efficient if one wants to use a 
Monte Carlo method to sum over all possible color 
states. The program will then extract through 
a standard unweighting (at the leading order in 
1/-/V?) a specific color flow from all possible color 
flows contributing to a given orthogonal color 
state. This color flow is then suitable as an initial 
condition for the shower evolution. The advan- 
tage of this approach is twofold: first the number 
of dual amplitudes contributing to the amplitude 
for external states with assigned colors is substan- 
tially smaller than the total one and second dual 
amplitudes are required only for accepted events 
which, in general, are a small fraction of the gen- 
erated ones. Further details can be found in [ |l|. 

3.1. Results 

As an example of our technique, we present 
here results for the following two parton-level pro- 
cesses: 

9 9 -> 8 g 

qq - 8g. (7) 

For comparisons, we also computed the above 
reactions in the simple approximation first sug- 
gested by Kunszt and Stirling [ |U| . This approx- 
imation (hereafter referred to as SPHEL) consists 
in assuming that the average value of maximally 
helicity violating amplitudes [ |ll| is equal to the 
average value of all other non-zero amplitudes. 

The kinematic configuration and the cut values 
used in our numerical examples are as follows: 

\/s = 1500GeV, p Tl >60GeV, \rn\ < 2, 
ARij > 0.7. 

These values, and the choice of a fixed strong cou- 
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Figure 1. Differential distributions for the 
minimum gluon transverse momentum. Ex- 
act result vs. SPHEL. 



Figure 2. Differential distributions for the 
maximun gluon transverse momentum. Exact 
result vs. SPHEL. 



pling a s = 0.12, only serve for illustrative pur- 
poses. 

In Fig. 1, we show the distribution of the min- 
imum gluon transverse momentum for both pro- 
cesses. In Fig 2 we plot the distributions for the 
maximum gluon transverse momentum. 

4. Conclusions 

We have reviewed the ALPHA algorithm to eval- 
uate the exact, tree-level matrix elements in the 
context of multi-parton processes in QCD. The 
algorithm is suitable (easy to implement) for the 
authomatic calculation of the scattering ampli- 
tudes and it allows the inclusion of mass effects in 
a straightforward manner. The CPU cost of the 
algorithm grows like K n (K constant, n number 
of external particles) as opposed to the factorial 
growth of the number of Feynman graphs. 

This technique has been tested for processes 
such as gg — > n gluons and qq — > n gluons, with 
n up to 9. We discussed how the summation over 



colour configurations allows the construction of 
parton-level event generators suitable to interfac- 
ing with a parton-shower evolution including the 
effects of colour-coherence. This will eventually 
lead to a fully exclusive, hadron-levcl description 
of multi-jet final states, accurately incorporating 
the dynamics of large jet-jet separation angles. 

REFERENCES 

1. F. Caravaglios, M. L. Mangano, M. Moretti 
and R. Pittau, Nucl. Phys. B539 (1999) 215 

2. F. Caravaglios and M. Moretti, Phys. Lett. 
B358 (1995) 332 

3. F. Caravaglios and M. Moretti, Z. Phys. C74 
(1997) 291; 

G. Montagna, M. Moretti, O. Nicrosini and 

F. Piccinini, Nucl. Phys. B541 (1999) 31; 

G. Montagna, M. Moretti, O. Nicrosini, 
A. Pallavicini and F. Piccinini, Nucl. Phys. 
B547 (1999) 39 



6 



4. M. Moretti, Nucl. Phys. B484 (1997) 3; 

G. Montagna M. Moretti, O. Nicrosini and 
F. Piccinini, Eur. Phys. J. C2 (1998) 483; 
F. Gangemi, G. Montagna, M. Moretti, 
O. Nicrosini and F. Piccinini, Eur. Phys. J. 
C9 (1999) 31; 

F. Gangemi, G. Montagna, M. Moretti, 

0. Nicrosini and F. Piccinini, Nucl. Phys. 
B559 (1999) 3 

5. J. C. Collins, Nucl. Phys. B304 (1988) 794. 

6. I. G. Knowles, Nucl. Phys. B310 (1988) 571. 

7. G. Marchesini and B. R. Webber, Nucl. Phys. 
B310 (1988) 461. 

8. G. Marchesini, B. R. Webber, G. Abbiendi, 

1. G. Knowles, M. H. Seymour and L. Stanco, 
Comput. Phys. Commun. 67 (1992) 465. 

9. G. Marchesini, B. R. Webber, G. Abbiendi, 
I. G. Knowles, M . H. Seymour and L. Stanco, 
hep-ph/9607393 . 

10. Z. Kunszt and W.J. Stirling, Phys. Rev. D37 
(1988) 2439. 

11. S.J. Parke and T. Taylor, Phys. Rev. Lett. 56 
(1986) 2459. 



